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Abstract 

Using the quasiclassical Green's function in the Coulomb field, we analyze the probabilities of 
single and multiple e + e~-pair production at fixed impact parameter b between colliding ultrarel- 



ativistic heavy nuclei. We perform calculations in the Born approximation with respect to the 



O 

' parameter Zbcx, and exactly in Z^a, Za and Zb are the charge numbers of the corresponding 



nuclei. We also obtain the approximate formulas for the probabilities valid for Zj^ol, Zbol < 1. 
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I. INTRODUCTION 



The cross section of e + e~ pair production in ultrarelativistic heavy-ion collisions is very 
large, and this process can be a serious background for many experiments. Besides, it is also 
important for the problem of beam lifetime and luminosity of hadron colliders. It means 
that various corrections to the Born cross section for one-pair production, as well as the 
cross section for n-pair production (n > 1), are very important. Recently, the process was 
discussed in numerous papers , see reviews jl|, |2j, |3( . However, some important aspects of 
the problem have not been finally understood, and in the present paper we are going to 
elucidate them. 

For our purpose, it is convenient to consider a collision of the nuclei A and B with the 
corresponding charge numbers Za and Zb in the rest frame of the nucleus A. The nucleus 
B is assumed to move in the positive direction of the z axis having the Lorentz factor 7. 
For 7 3> 1, it is possible to treat the nuclei as sources of the external field, and calculate the 
probability of n-pair production P n (b) in collision of two nuclei at a fixed impact parameter b. 
The corresponding cross section a n is obtained by the integration over the impact parameter, 

°n = J d 2 bP n (b). (1) 

Average number of the produced pairs at a given b reads: 

W(b) = J2 nP n(b)- (2) 

n=l 

The function W(b) defines the number- weighted cross section 

/oo 
d 2 bW{b) = J2 nCT n- (3) 
n=l 



The closed expression for o T was obtained in Refs. m , though the correct meaning of 
this expression was recognized later in Ref. p|- 

The cross section <jt can be presented in the form: 



(JT — ®T ' °T ' ®T 1 V±J 

where 0% is the Born cross section, i.e., the cross section calculated in the lowest-order per- 
turbation theory with respect to the parameters Z^ B a ( 0% oc (Z B ot) 2 (Z AOt) 2 , a = e 2 is the 



2 



fine-structure constant, e is the electron charge, K = c = 1), cr^ is the Coulomb corrections 
with respect to one of the nuclei (containing the terms proportional to {Zbci) 2 {Z A a) 2n or 
(ZB<y) 2n {Z A a) 2 , n ^ 2), and a^ c is the Coulomb corrections with respect to both nuclei 
(containing the terms proportional to (ZB<y) n {Z A a) 1 with n,l > 2). The cross section Oj, 
coincides with the Born cross section of one pair production, which was calculated many 
years ago in Refs. 0, 0] . 

The expression for W(b) derived in Refs. jj, 0, |(| requires regularization. The correct 



regularization was made in Refs. 



10 



111 ], where the expressions for er^ and <7f. were 



obtained in the leading logarithmic approximation: 
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[Z A a) 2 (Z B a) 2 , L 



~^^- 2 L 2 [f(Z B a) + f(Z A a)} 
97T m l 
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9tt m 
In 7 



2 Lf(Z B a)f(Z A a) 



f(x) = Re[ip(l + iZ A a)+C) 



(5) 



where m is the electron mass, ip{x) = T'(x)/T(x), and C = 0.577. . . is the Euler constant. 
The expression for coincides with that obtained in Ref. 12j by means of the Weizsacker- 
Williams approximation. The accuracy of the expression (@J) with cr^ and a^P given in 
(0) and a Bom from Refs. j^, 0] is determined by the relative order of the omitted terms 
~ [Z a ,boi) 2 I L 2 . This accuracy is better than 0.4% for the RHIC and LHC colliders. In the 



recent papers 



13 



14 1 , the Coulomb corrections were calculated numerically for a few values 

nn 



of 7. We emphasize that the accuracy of the results in Refs. 13j, [l4j is the same as in (J^J). 
The uncertainty is related to the contribution of the region, where the energies of electron 
and positron are of the order of electron mass in the rest frame of one of the nuclei. 

In Refs. H H Q it was Canned that the factorization of the m nltip, pan pro- 
duction probability is valid with a good accuracy, resulting in the Poisson distribution for 
multiplicity: 

W n (h\ 

(6) 



nl 

The factor exp(— W) is nothing but the vacuum-to-vacuum transition probability Pq = 
1 — Y^m=i Pn- Strictly speaking, the factorization does not take place due to interference 
between the diagrams corresponding to the permutation of the electron (or positron) lines 
(see, e.g., 3]). Nevertheless, one can show that this interference gives the contribution which 
contains at least one power of L less than that of the amplitude squared. Therefore, in the 



3 



leading logarithmic approximation one can use the expression (jHJ). Thus, to obtain a n it is 



sufficient to know t 
In Refs. |lS, 20, 



re function W(b) 



21 



22, 



231 ] . the function Wo(b) ( the Born approximation for W(b)) was 



calculated numerically for mb < 1 and a few particular values of 7. The correct dependence 
of the function Wo(b) on b at mb 3> 1 was obtained analytically in Ref. 2M by two different 
methods. Both methods give the following result: 

28 C 



Wo{b) ^(mbf 
in the region 1 <C mb < */7, and 

W (b) 28 



[2 In 7 — 3 In (mb)] In (mb) 



c 



mb) 



(7) 



(8) 



9tt 2 (mb) 2 

in the region yfiy < mb <C 7. Note that the function Wo(b) given by Eqs. (|IJ) and (JBJ) is the 
continuous function at mb = 7 together with its first derivative. Certainly, the integration 
of Wo(b), Eqs. (|7j) and (jH)) over b gives the main term (oc L 3 ) in cr^. In the recent paper 
. the ansatz for W (b) has been suggested, that has quite different dependence of W (b) 
on 7 and b for 1 <C mb <C yfy. In the present paper, we confirm the result (JJJ) once more 
and unambiguously disprove the ansatz suggested in Ref. 2^|. 

The cross section o\ of one-pair production can be represented as follows: 



o T + a unit = j d 2 bW(b) - j d 2 bW(b) (1 - e~ w ®) . 



(9) 



Thus, the difference between <j\ and <jt is due to the unitarity correction <T un i t . The main 
contribution to the term Ot comes from b 1/m. It was shown in Ref. 2M that the main 
contribution to the second term, cr unit , as well as the main contribution to the cross sections 
for the n-pair production (n ^ 2), comes from b ~ 1/m. As shown in Ref. 2J|, in this 
region, the function W(b) has the form 



W(b) = C,LT(mb) 



(10) 



where the function T(mb) depends on the parameters Z^a and Z^a and is independent of 
7. Let us represent the function J-(x) as 



F(x) = Fq(x) + Ta{x) + Tb{x) + Tab(x) 



(11) 



where T§(x) is independent of Za and Zb (Born term), Ta(%) contains terms oc 
(Z A<y) n>0 (Z bo) (Coulomb corrections with respect to the nucleus A), JF B (x) contains terms 



oc (Z a®) c \Z B<y) n>0 (Coulomb corrections with respect to the nucleus B), and Fab(x) con- 
tains terms oc (Z A<y) n>c \Z b«) i>0 (Coulomb corrections with respect to both nuclei). 

In the present paper, we calculate the function J-{x) for Z^a <C 1 , Z^oi < 1, and x < 1. 
In this limit, one can neglect in Eq. the terms J-'six) and Tab{x). Though Z^a <C 1, 
we can not expand exponent in ((BJ) if £L ~ 1. Our method is based on the use of the 
quasiclassical Green's function of the Dirac equation in the Coulomb field. 



II. GENERAL DISCUSSION 

In the leading in Zga order, the matrix element of e + e~ pair production, M, has the 
form 

M = -e J dtdr exp[-i(e p + e q )t]H p _(r) A{t, r)^_ p+ (r), (12) 

where *4. M (t, r) is the four- vector potential of the moving nucleus B, ty p _ and ^- p+ are 
the positive- and negative-energy solutions of the Dirac equation in the Coulomb field of 
the nucleus A; p_ = (e p , p), p+ = (e q , q) are the four- momenta of electron and positron, 
respectively. 

Then we use the Fourier transform, A%, of the vector potential A^{t, r), 

< = -kiTWTW^^ 6 {lk ° ~ lPk * ] u " ' (13) 

where = (7, 0, 0, 7/3) is the four-velocity of the nucleus B, and b is the impact parameter. 
Taking the integrals over t, k°, and k z , we obtain 

An 7 ry C dk p —ik ± -b r 

M = I J^ WheJW* J * ^i^^^VMW-.W, (14) 

where E = e p + e q , r = (p, z). 

At the calculation of the probabilities integrated over the angles of the final particles, 
it is convenient to exploit the Green's functions of the Dirac equation in the external field. 
Using the relations (see, e.g., j3]) 

\2 



[ dQ q ^_ p+ (r 2 )^_ p+ (r!) = -i^-8G {r 2 , r x \ - e q ) , 

[ dQ p %_(r 1 )%_(r 2 )=i^-SG(r 1 ,r 2 \e p ), (15) 
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where SG (r, r'\e) is the discontinuity of the Green's function on the cut, and the summation 
is performed over the spin states, we obtain for the total probability: 



W(b) = J2\M\ 



, dpdq 

(2^)6 



2Zboc\ f de q de p dkx±dk2± 

~^F) J (27T) 4 



exp [i(k 



i± 



k 



2±j 



[k 2 1± + ( W) 2 ] m ± + {Eh(3f 



ik 2 ± • p 2 - ik 1± • pi + - -I. : 



x J dr\dr 2 exp 
x Sp [u8G(r 2 , r x \ - e q ) u 5G(r x , r 2 \e p )\ . 



.E 



(16) 



Using the gauge invariance and the condition 7 ^> 1, it is possible to make in Eq. (fT6j) the 
following replacement 



7 

Sp [u5G{r 2 ,ri\ -e q )u 8G{r x ,r 2 \e p )\ -> — Sp 



A; 2 ±5G(r 2 ,r 1 | - e q ) k 1± 5G(r 1 ,r 2 \e p ) 



(17) 



In the leading logarithmic approximation, the main contribution to the probability W(b) 
comes from the region e± ^> m, where the quasiclassical approximation is applicable. Be- 



sides, it is convenient to perform t 



of the squared Dirac equation 
we obtain 



25 



re calculations in terms of the Green's function D(r, r'\e) 
Using the transformations similar to those in Ref. 



W(b) = 4 (Z B af 



2 f de q de p dki±dk 2 ± exp [i(k\± — k 2 ±) • b] 

E* (2tt) 4 W± + WW] IK± + ( E hW 

E 

ik 2 ± • p 2 - iki± • p x + l ~p{ z 2 ~ z i) 



x J dr\dr 2 exp 

x Sp I [[-2ifc 2 ± • V 2 + k 2 h±] D(r 2 , n\ - e q 
x [[-2ifei_L ■ Vi - fcifeij.] D(n, r 2 \e p 



(18) 

Here k\ = (E, k-\±,E), k 2 = (E, k 2 ±, E). In the quasiclassical approximation, the function 
D has the form |25j ] 



D(r 2 ,ri\e) 



8n 2 rir 2 



dq exp 



nr(q + ff 



2r 1 r 2 



V Q 2 J 



X 



Xr 



l + o <*-(q + f) 

2r x r 2 



K 



A 



K 



a = 7°7 , / 



[[n x r 2 ] x r] 



r = n-r 2 , (19) 



where q is a two-dimensional vector lying in the plane perpendicular to r. The explicit form 
()19|) of the quasiclassical Green's function is very convenient for analytical investigation of 
high-energy processes in the Coulomb field. 



III. ANALYTICAL RESULTS 



For mb < 1, the main contribution to the integrals in Eq. (fTBj) is given by the re- 
gion of small angles between vectors r±, —r 2 , and z-axis. Using these conditions and the 
quasiclassical Green's function (|19|). we obtain the following representation for T{mb) = 
jF (m6) + jF^mfr) (details of the calculation are presented in Appendix) 

\R + xQ\^ 2lZAa ~ 



^(mb) 



vr 4 (Z A a) 



o 



R-xQ 

x UV^{x-x)l3-Q[Kl{Q)/Q 2 -K 1 {Q)K 1 {Q)/QQ) + \K (Q) - K {Q) 
+Axx(3 2 K 2 (Q)/Q 2 + (Q 2 - Axx (J3 ■ Q) 2 / (3 2 ) [#i(Q)/Q - #i(Q)/q] * 
Q 2 = Q 2 + f3 2 , R= Vxlcp + mb, x = l-x, 

(20) 

where K n (x) is a modified Bessel function of the third kind. The form (|2Tlj) is suitable for 
investigation of the asymptotics of !F{mb). For numerical evaluation, it is convenient to pass 
from the integration over the angle of the vector Q to the integration over the parameter 
v using the following identities 

1 

COS0 

cos 20 



[d<j> 


1 - 


(r 


+ a 


COS0\ 


J 2tt 






- b 


COS0 / 



v sinh txv 1 
hm 

7T (5^0 /q 



dv 



V 1 5 V 1 s \v 



In 



l+VT- 



2 \i+VT- 



v = 1 — v , s = av — bv . 
Making the substitution v = u/(u + where 

u = l-u, f = 



(21) 



R 2 + x 2 Q 2 
R 2 + x 2 Q 2 ' 



7 



and taking into account the symmetry of the integrand with respect to the substitution 
u — > u, x — ► x, we obtain 



^(mb) 



n A Z A a 



x ( In 



+ VI 



C a r Ann f d2 P f 1/2 du \y w i-w 
f ctx / aQJQ / — - / — r cos [Z^am (m/mJJ 

o Jo J P Jo uu 



l-2xx) QK 1 (Q)/Q-K 1 (Q) 



t 2 



+ 



K (Q) - K {Q) + 4xxK 2 (Q)p 2 /& 



xx 



(3 R 

I3R 



-Ayfxx (x — x) 



t 2 



(3 R 
R 



RQx 



-\ 2 



QK^/Q-K^Q) 
RQx 



t 2 



9 



QKf(Q)/Q 2 -K 1 (Q)K 1 (Q)Q 



Q 2 = Q 2 + f3 2 , g = max( J R 2 , Q 2 x 2 ) , 
t = 2QR (xu - xu) , s = R 2 + Q 2 (x 2 u + x 2 u) 



(22) 



Let us consider the asymptotics of Eq. (J2(J|) . For m& ^> 1, there are two regions of 
integration over (3, giving the leading logarithmic contribution to J- {rob): 1 <C \\fxx T /3 + 
mb\ <C mb and 1 <C (3 <C mfo. These regions give equal contributions, and the final result 
reads 



T(mb) 



56 



In (mb) . 



(23) 



97r 2 (mb) 2 

Thus, the leading logarithmic contribution is given by the Born term To(mb). This asymp- 
totics agrees with Eq. (jJJ) under the condition ln(mfe) <C L. 

The main contribution to T^mb) comes from the region \\Jx!cf3 + mb\ ~ 1 and has the 
form 



^A(mb) 



28 



(24) 



97r 2 (m6) 2 ' 

where the function f(x) is defined in Eq. (jSJ. Again, this asymptotics is valid under 
the condition \n(mb) <C L. Similar to the derivation of Eq. (JZj), see Ref. [2^, based on 
the equivalent photon approximation, it is possible to obtain for the Coulomb corrections 
to W(b) with respect to the nucleus A, Wa(J>), the expression valid in the wider region 
ln(m&) < L (but still 1 Cm6<7). We have 



W A (b) 



-|^W /( ^ a)ln (^) 



(25) 



Eq. pij) evidently agrees with Eq. (J2SJ). 

Let us consider the asymptotics at small impact parameters. For mb <C 1, the leading 
logarithmic contribution comes from the region mb <C (3 ~ Q -C 1. Taking the integrals 
over this region, we obtain 

T(mb) — I In f ) Re + iZ A a) + C - (Z A a) 2 

37T 2 {Z A a) 2 \mbj (26) 

+iZ A a(l + (Z A a) 2 )^'(l + iZ A a)] . 

This asymptotics is obtained for zero nuclear radius R n . In order to obtain W(b) for the 
extended nuclei, it is sufficient, within the logarithmic accuracy, to make the substitution 
In (mb) — > In (mb + mR n ) in the asymptotics (|26|). For 6 ^> i? n , the finite-nuclear-size 
correction to is negligible. 



IV. NUMERICAL RESULTS 

Using Eq. ()18jl . we performed the tabulation of the function T(mb) for a few values of 
Za- The corresponding results are presented on the left plot of Fig. ^ and in the Table HJ 
We remind that these results are obtained in the Born approximation with respect to the 
nucleus B. For most experiments Z A = Zb, and it is necessary to know the function T(mb) 
beyond the Born approximation with respect to the nucleus B. If we assume that the term 
Tab m Eq. (Jll|) is numerically small, then we can approximate the function T as To + IT a 
in this case. This function is shown on the right plot of Fig. ^ It is seen that the Coulomb 
corrections in the region mb < 1 are very important for the experimentally interesting case 
Za = Z B = 79. The assumption of smallness of the contributio n Tar i s supported by the 



comparison of our results for W(b) with those obtained in Refs. 14L \27\ for Z A = Zb = 79 
and 7 = 2 x 10 4 (j c . m . = 100). 

As we already pointed out, Eq. (jlOj) has logarithmic accuracy which can be sufficient for 
very large 7. In order to go beyond the logarithmic accuracy, we represent W(b) in the form 



W(b) =([L- G(mb)} T(mb) , (27) 

where G(mb) is some function of mb and, generally speaking, of the parameters Z A a and 
Z B oi. The asymptotics of G(mb) at 1 <C mb ^/j is known, see Eqs. flj|) and (J23j). 
However, the calculation of the function G(mb) at mb < 1 is rather complicated problem. 



9 



nstead, we use the results of numerical calculations, performed for a few values of 7 in Refs. 



19 



271 ] in the Born approximation. We have found that the form 



G{mb) = -ln(m& + 1) + 1.9, 



(28) 



19 



27 



111 



provides good agreement of Eq. (}2Tj) with the numerical results of Refs 
the wide region of mb, see Fig. |2j The form (}28|) of G(mb) is obtained by fitting the Born 
results and thus is independent of Za,b- It provides the correct asymptotics of Wo(b), Eq. 
(JZJ). It turns out that the formula (f2*Tj) with G(mb) from Eq. (|2*%|) has also high accuracy 
for Z^a, Zbol < 1 in the region mb < 1 where the Coulomb corrections are large. We have 
checked this fact by comparing our results with those of Ref. (27 1 obtained numerically for 
Z a = Z r = 79, see Fig. El Note that the tabulation of W(b) and P/v(&), performed in 

nn nn n 

Refs. [121, HJ, [IjJ, |27J, |2SJ] for a few values of 7, required the evaluation of nine-fold integral 
and, therefore, was very laborious. The calculation of the function T from Eq. (|2*2*|) is 
essentially simpler. Besides, since this function is independent of 7, one can easily obtain 
predictions for W{b) at any 7 ^> 1 using Eqs. (j2*7|) and 



V. CONCLUSION 



In the present paper, we have found in the leading logarithmic approximation the simple 
representation for the function W(b) for mb < 1, Z^a <C 1, and arbitrary Z^a. Using the 
results of numerical calculation of W(b) performed for a few values of 7 and Z^.b, we have 
obtained the approximate formula for W(b) valid in a wide region of parameters: mb < */7, 
Zaol < 1, Zb« < 1, 7 ^> 1. We estimate the accuracy of this formula to be a few percent. 
The results obtained clearly demonstrate the dependence of W{b), as well as P n (b), on the 
relativistic factor 7 and the parameters Za,b&- 

This work was supported in part by RFBR grant No. 05-02-16079 and by the grant for 
young scientists of SB RAS (R.N.L.). 



APPENDIX A: CALCULATION OF THE INTEGRALS 

In this Appendix, we present some details of derivation of Eq. (}20|) from Eq. (|T8|) . The 
main contribution to the integrals comes from the region of small angles between vectors r%, 
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TABLE I: The function J-(x), Eq. (|22[). calculated in the Born approximation (Zj^a — ► 0) and 
exactly in the parameter Zj^a for Au, Pb, and U. 

— V2, and z-axis. Using this fact, we take the integrals over the angles of r% and V2, make the 
substitution —>■ Erip, and change the variables e p = Ex, e q = Ex = E(l — x). Taking 
the integral over E in the logarithmic approximation at 7 3> 1 and mb < 1, we obtain 
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FIG. 1: The function F(x), Eq. $V$, for Z A = 92 (dash-dotted line) , Z A = 79 (dotted line), 
Za = 47 (dashed line) , and the Born approximation (solid line) . Left plot corresponds to the Born 
approximation in Zb<x, Eq. (|22|) . Right plot shows the results obtained from Eqs. and (|22l) 

for Zb = Za with the term Tab (%) omitted. 



3.5 




dW(b)= { -^t In 7 



(2vr 
x exp 



dki± dk 2 ± 

~P P 

"-i_l ^2± 



dxxx — - — - / dQ dq 
ri r 2 



|g + QI 
|g-QI 



- 2 2 ^(r 1 + r 2 )Q 21 

—m (ri + r 2 ) — i A ■ p xa^ri/c-^ + r 2 k 2± ) + ■ 



\i 2 {ri + r 2 ) - i A • /3 - - 
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r 2 
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rir 2 
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r 2 fc 2 ± 



2xa;rir 2 2r 2 2r x 

A = fc 1± -fc 2± , P = q/2 + (x-x)Q/2-b. 
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FIG. 2: The probability of one-pair production -Pi (6) corresponding to the function W(b) from Eq. 
(|27|). 7 = 2 x 10 , Za = -^b = 79. Dashed line: the function T is taken in the Born approximation, 
T = J-q; solid line: the Coulomb corrections are taken into account, J- = J-q + 2J 7 a- Dots show 
the corresponding results of numerical calculations from Ref. |27| . 



The integration over two-dimensional vectors ki± and k 2 ± can be easily performed. The 
result reads 

'' '■• ' '' i{r l + r 2 )Q 2 



dW(b) 



3 J In 7 I dx—^—^- / dQdqexp 



x 



+ 



(27T) 

\q + Q\ 
\q-Q\ 



2iZ A a 



n r 2 



2 [x — x) 



P Q , i 



+ 



4(/3-Q) 2 m^rx + ra) 2 
xa; + - 



/3 4 

£j = exp[i(3 2 /2xxr 



2(3 Z 



f3 2 ' 2xx 
{S 1 - 1) [E 2 - 1 



-m (ri + r 2 ) + 
2 2rir 2 

(2^ -E x - £ 2 ) - 4£ 1( ? 2 



(A2) 
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Taking the integrals over rip and passing from the variable q to f3 = q/2 + (x — x) Q/2 — b, 
we obtain Eq. (QUI) . 
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